Momentum distributions in medium and heavy exotic nuclei 
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I. INTRODUCTION 

The study of nuclei close to the nuclear drip line 
and even beyond it has been greatly extended in recent 
years. This increased interest is based on new phenom- 
ena that already have been observed or that are predicted 
to occur in these nuclei. Since the first experiments 
[11, S S 0, S @1 , it has been found from analyses of total 
interaction cross sections that weakly-bound neutron-rich 
nuclei, e.g. ^'^He, ^^Li, ^''Be, ^^'^^B, have increased sizes 
that deviate substantially from the R ~ A^^^ rule. It was 
realized (e.g. Refs. 7, 8, 9]) that such a new phenomenon 
is due to the weak binding of the last few nucleons which 
form a diffuse nuclear cloud (the "nuclear halo") due to 
quantum-mechanical penetration. Another effect is that 
the nucleons can form a " neutron skin" Q when the neu- 
trons are on average less bound than the protons. The 
origin of the skin lies in the large difference of the Fermi 
energy levels of protons and neutrons so that the neutron 
wave function extends beyond the effectively more bound 
proton wave functions Q. 

The experiments on scattering of radioactive nuclear 
beams on proton target at various incident energies (e-g., 
less than 100 MeV/nucleon for He isotopes p^. [ill. fisl. 
H, M, M, PlQEI and 700 MeV/nucleon for He and 
Li isotopes |20l. l2ll 122. [23l. [2^ ) have allowed one to study 
the charge and matter distributions of these nuclei using 




As known, the most accurate determination of the 
charge distributions in nuclei can be obtained from 
electron-nucleus scattering. For the case of the unstable 
exotic nuclei the corresponding charge distributions are 
planned to be studied by colliding electrons with these 
nuclei in storage rings (see, e.g. the G SI p hysics pro- 
gram [39|] and the plan of RIKEN ^ |41[). A num- 
ber of interesting issues can be analyzed by the electron 
experiments. One of them is to study how the charge 



distribution evolves with increasing neutron number at 
fixed proton number, or to what extent the neutron halo 
or skin may trigger sizable changes of the charge root- 
mean-squared (rms) radius and the diffuseness in the pe- 
ripheral region of the charge distribution. In Ref. ^45| we 
studied charge form factors of light exotic nuclei (^'^He, 
"Li, "Be, "d9B) using various theoretical predictions 
of their charge densities. Among the latter we used those 
from Tanihata et al. @ for the He isotopes, those from 
the cluster-orbital shell-model approximation (COSMA) 
for the He [2^ and Li [l^ isotopes, those from the mi- 
croscopic large-scale shell-model (LSSM) method (for He 
[i^ l and Li and those from Suzuki et al. [45*1 for 

"'^^Be and ^'^'^^B nuclei. In Ref. 46.] our calculations of 
the charge form factors of exotic nuclei were extended 
from light (He, Li) to medium and heavy nuclei (Ni, Kr, 
and Sn) . For the He and Li isotopes the proton and neu- 
tron densities obtained in the LSSM method have been 
used, while for Ni, Kr, and Sn isotopes the densities have 
been obtained in the deformed self-consistent mean-field 
Skyrme Hartree-Fock (HFUBCS method [13, IH, |4|. 
In contrast to the work [Z^l, in [i^ we calculated the 
charge form factors not only within the PWBA, but also 
in DWBA by the numerical solution of the Dirac equa- 
tion [5^, mi, [13 for electron scattering in the Coulomb 
potential of the charge density of a given nucleus. The 
role of the neutrons has been taken into account. A de- 
tailed study of the charge radii and neutron skin in Ni, 
Kr, and Sn nuclei, as well as of the formation of the pro- 
ton skin, has been performed within the same method in 

Another important characteristic of the nuclear ground 
state is the nucleon momentum distribution (NMD) n{k). 
The scaling analyses of inclusive electron scattering from 
a lar ge v ariety of nucleifsee, e^ [s^lssj for the y-scaling 
and j56l . IstI . IssI . IsqI . [60 . [6ll . WX l63l | for i/)'-scaling and 
superscaling) showed the evidence for the existence of 
high-momentum components of NMD at momenta k > 2 
fm^^. It has been shown [b^, [61|, [13, [l^ that it is due 
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to the presence of nucleon-nucleon (NN) correlations in 
nuclei (for a review, see e.g. [13]). It has been pointed 
out that this specific feature of n{k)/A is similar for all 
nuclei, and that it is a physical reason for the scaling and 
superscaling phenomena in nuclei. As known 64. i65l|. the 
mean-field approximation (MFA) is unable to describe si- 
multaneously the two important characteristics of the nu- 
clear ground state, the density and momentum distribu- 
tion. Therefore, a consistent analysis of the effects of the 
NN correlations on both quantities is required using theo- 
retical methods beyond the MFA in the description of rel- 
evant phenomena, e.g. the scaling ones. Particular atten- 
tion has been paid to the NMD in a given single-particle 
state analyzing the (e, e'p) reactions in nuclei (see, e.g. 
the review in '66||, the works and ref- 

erences therein). The self-consistent density-dependent 
HF (DDHF) approximation has been applied in 48] to 
calculate NMD in spherical and deformed Nd isotopes, 
studying the effects of deformation, as well as those of 
pairing and of dynamical short-range NN correlations. 

It is of importance to study the NMD not only in 
stable, but also in exotic nuclei. It is known (see, e.g. 
Ref. [72|) that in the reactions with exotic nuclei the mo- 
mentum distribution of a core fragment of the projectile 
reflects the momentum distribution of the valence nucle- 
ons in the projectile near the surface [73|. Many exper- 
imental (e.g, and theoretical [tHI works have been 
carried out to study the momentum distribution from 
the breakup of the projectile. For instance, in Ref. [t^ 
the momentum distribution of relative motion between 
two nucleons has been calculated for both ^He and ^Li 
two-neutron halo nuclei. The obtained results in the case 
of realistic NN interaction show two interesting predic- 
tions: i) S'-wave dominance in the NMD of ^He and ii) 
the ^Li momentum distribution is very similar to that 
of the deuteron. In fll\ the neutron and proton mo- 
mentum distributions in some stable nuclei (^^C, ^^O, 
''"Ca, ^^Fe, and ^°*Pb) were calculated along with those 
of hght neutron-rich isotopes of Li, Be, B, and C using 
the natural-orbital representation (NOR) on the basis of 
the empirical data for n(fc) in ^He and, especially, those 
for the high-momentum components of the latter. 

The main aim of our work here is to calculate the NMD 
for the same isotopic chains of neutron-rich nuclei (Ni, 
Kr, and Sn) for which we had studied charge densities, 
radii, form factors, halo, and skin in our previous works 
and [11]. The mean- field contributions to n{k) in 
these nuclei are calculated within the same self-consistent 
approach applied there in which the one-body energy 
density functional is obtained starting from a two-body 
density-dependent Skyrme interaction and a pairing in- 
teraction that is treated in the BCS limit. The HF equa- 
tions are solved for the {N, Z) nucleus using a deformed 
harmonic-oscillator basis in cylindrical coordinates with 
oscillator lengths used as variational parameters. The 
BCS equations are solved at each HF iteration and the 
occupation numbers are used to construct the density- 
dependent mean field for the next HF iteration. We refer 



to this mean field approach as DDHF-I-BCS. The remain- 
ing effects of the NN interactions, to which we refer as 
NN correlations, are considered in two of the correlations 
approaches, namely, in the approach (see [6l|, [H, ^7^) 
using the light-front dynamics (LFD) method (e.g., [79j ) 
and in that [sO] based on the local density approxima- 
tion (LDA). Several questions are investigated, such as 
the sensitivity of n{k) to all details of the calculations, 
e.g.: i) to different types of Skyrme forces; ii) to the 
pairing correlation effects; iii) to the effects of nuclear 
deformation; iv) to the strength of the NN correlations 
included in the LFD and LDA approaches (respectively, 
to the values of the correlations strength parameters /3 
and 7). A special attention is paid to the isotopic and 
isotonic sensitivity of the proton and neutron momentum 
distributions. The results for n{k) in the exotic nuclei are 
compared with that in nuclear matter (NM). 

The paper is organized in the following way. Section HIl 
contains the formalism of the deformed Skyrme HF-I-BCS 
method and the approaches based on the LFD and LDA 
methods that provide the model nucleon momentum dis- 
tributions. The numerical results and discussions are pre- 
sented in Sec. IIIII Finally, we draw the main conclusions 
of this study in Sec. IIVI 



II. THEORETICAL FRAMEWORK 
A. Deformed Skyrme HF+BCS formalism 

Some of the results discussed in the next Section have 
been obtained from self-consistent deformed Hartree- 
Fock calculations with density dependent Skyrme inter- 
action [49] and pairing correlations. Pairing between like 
nucleons has been included by solving the BCS equations 
at each iteration either with a fixed pairing gap param- 
eter (determined from the odd-even experimental mass 
differences) or with a fixed pairing strength parameter. 
We consider in this paper the Skyrme force SLy4 [8l| 
which gives an appropriate description of bulk properties 
of spherical and deformed nuclei. 

Assuming time reversal invariance, the single-particle 
Hartree-Fock solutions for axially symmetric deformed 
nuclei are characterized by the eigenvalue fli of the third 
component of the total angular momentum on the sym- 
metry axis and by the parity tt^. The state i can be 
written as 

+ <f>r(r^,z)e^A^^;,_(a)], (1) 

where Xqi{Q)j X±{'^) ^-^e isospin and spin functions, A^ — 
i7i ± 1/2 > 0. r±, z,ip are the cylindrical coordinates of 
r. 

The wave functions are expanded into the eigen- 
functions (f>a of an axially symmetric deformed harmonic- 
oscillator potential in cylindrical coordinates. We use 12 
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major shells in this expansion, 

(r, CT, q) = Xq. (q) ^) 



with a = {n±, Uz, A, S} and 



6a (r,cr) = V„^(rx)V'„,(z)--^Xi:(c^), 



in terms of Hermite and Laguerre polynomials 



(2) 



(3) 



(4) 



71 _L 



K +A)! 



/3^V2 7^A/2e-"/2<(,7) (5) 



with 



(6) 



From the expansion ^ we may conveniently express 
the single-particle Hartree-Fock wave functions in mo- 
mentum space, which we denote as ct, q). They are 
given by [48j 



(7) 



with 



4>a{k, a) = 



(27r)3/2 



dre-''^-^>„(r,(T) (8) 



normalized to one. 

The spin-independent proton, neutron and total den- 
sities are given by 



p{r) = p{r^,z) = ^ 2ufp,( 



(9) 



in terms of the occupation probabilities vf resulting from 
the BCS equations and the single-particle densities pi 

p,{r) = p,(rx, z) = |<i>+(ri, z)\' + \^;{r^, z)\' , (10) 

with 
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2tt 



E '5s,±i/2'5a,at C<(r±)V'n.(4ll) 



Similarly, we define in momentum space the proton, neu- 
tron and total momentum distributions by 

n{k) ^ n{k^,kz) ^Y.2v^n,{k±,k,) , (12) 



where k±,kz are the cylindrical coordinates of k. The 
single-particle momentum distributions ni{k) are given 
by 



n,ik) = n,{k^,kz) = |$+(fci,A:,)|2 + |$-(A:i,fc,)|2 
where 

^Hk±,kz) = 



/27r 



(13) 



(14) 



E '^s.±i/2'5a,atC;(-z)^V^^ (A;,)V'„, (fc«). 



fc2 _|_ ^2 

fc^ = kz/Pz, kri = k^/ f3j_ ~ "2 , 



(15) 



and = A+nz + 2n± is the major shell quantum number 
of the basis state a. 

The multipole decomposition of the momentum distri- 
bution can be written as 

n{k) = J2nx{k)Px{cos9k) 

A 

= no{k)+n2{k)P2{cos9k) + ... (16) 
with multipole components A 



nx{k) = 



2A + 1 



Fx (cos 0k) 



X Ti{k cos 9 k, k sin 9 k)d{cos 9k). (17) 



B. Methods going beyond the mean- field 
approximation 

It is well known (e.g. [g^l) that the methods within 
the MFA (e.g. shell model, Hartree-Fock and others) can 
describe the nucleon momentum distribution n{k) only 
for momentum values fc <1.5 fm^^ and are unable to ex- 
plain n{k) for larger k. The high-momentum components 
of n{k) {k >1.5 fm~^) are due to the specific forces be- 
tween the nucleons near the nuclear core (rc ~ 0.4 fm) 
that are the reasons for the short-range and tensor NN 
correlations. The differences between the values of n{k) 
for large k obtained within the correlation methods (e.g., 
exp(S')-method [83 |. the Jastrow correlation method and 
others, for a review see e.g. ^^) reach orders of mag- 
nitude. In this Subsection we consider the effects of NN 
correlations included in two correlation methods on the 
high-momentum contributions to the nucleon momentum 
distribution. 



1. Theoretical approach based on the light-front dynamics 
method 



Here we derive the NMD obtained within the LFD ap- 
proach developed in Refs. [b!], [T^I . The latter is based 
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on the nucleon momentum distribution in the deuteron 
from the Hght-front dynamics method (e.g. [zl])- Using 
the natural-orbital representation of the one-body den- 
sity matrix 'sS], n{k) is written as a sum of the contribu- 
tions from the hole-states [h'^{k)] (states up to the Fermi 
level (F.L.)) and the particle-states [n-^(fc)] (see also (6ll |) 
for protons (Z) and neutrons (N): 



(18) 



the well- 



where 



Concerning the particle-state contribution In" 
in (UHl) and ^ we used in Refs. ^ ^ 
known facts that: (i) the high- momentum components 
of n(fc) caused by short-range and tensor correlations 
are almost completely determined by the contributions 
of the particle-state natural orbitals (e.g. [Ill), and (ii) 
the high-momentum tails of the momentum distributions 
per particle are approximately equal for all nuclei and are 
a rescaled version of the nucleon momentum distribution 
in the deuteron nd{k) [ssl-ls^. 



and 



C{k) 
Z{N) 



C{k) 
Z{N) 



F.L. 



^2(2j + l)A„,,|i?„i,(fc)p (19) 



nlj 



Y,2{2j + \)Kaj\RraAk)?- (20) 



F.L. 



In (Ull and m 



UAik) ~ C^ridik), 



(26) 



where C is a constant. These facts made it possible to 
assume in |6TI.[78| and later, using the modification of the 
approach in [63|, that n^^^j.^-^{k) is related to the high- 
momentum components of the nucleon momentum dis- 
tributions in the deuteron. Thus 'nFz(N) i^) from Eq. ([201) 
can be substituted by (up to a normalization factor): 



C(fc) 



ruN 



(21) 



niN being the nucleon mass, Xnij are the natural occu- 
pation numbers and R„ij{k) are the corresponding wave 
functions in fc-space of protons (neutrons) in states with 
quantum numbers nlj . The momentum distribution (11^ 
is normalized to unity: 



nz(N){k)dk = 1. 



(22) 



As shown in [78|, in the NOR the hole-state contri- 
bution n^(^j^-^{k) to the momentum distribution can be 
represented to a good approximation by the momentum 
distribution from the spherical mean-field methods, i.e. 
\nij are taken to be unity and Rnij{k) are replaced by 
the corresponding mean- field cigcnfunctions. In our work 
(fc) in (0 and (HH) by 



we substitute n^z{N) 



Z(N) 



{k) = 



C{k) . 



Z{N) 



nz(N){k), 



where nz(N){k) is expressed by the NMD n^l^^ (k) ob- 
tained within the DDHF formalism (see Sec. II. A.) 



nz{N){k) 



^(^)^Z(^) (fc) 
/dfc'C(fc')nf(^f^(fc') 



and the normalization is 



C{k)fiz(N){k)dk = Z{N). 



(24) 



(25) 



We note that rj-l-j-jy-) (fc), nFz(N)(k) and nz(N){k) include 
the function C(fc) that originates from the relativistic 
LFD approach. 



"z(Ar)(^) = P[n2{k) + n^ik)] , 



(27) 



where /3 is a parameter, and n2{k) and n^{k) are ex- 
pressed by angle-averaged functions fl^ as: 



n2{k) = C{k)fl{k) 



(28) 



and 



n,{k)^C{k){l-z^)fl{k). 



(29) 



In (|^ z = cos{k,n), n being a unit vector along the 
three-vector (w) component of the four-vector to which 
determines the position of the light-front surface [79j . 
The functions f2{k) and /5(fc) are two of the six scalar 
(23) functions fi-e{k^,n- k) which are the components of the 
deuteron total wave function ^{k,n). It was shown [t^ 
that /s largely exceeds other /-components for k > 
2.0-2.5 fm~^ and is the main contribution to the high- 
momentum component of nd{k), incorporating the main 
part of the short-range properties of the NN interaction. 
It was shown in 63], however, that not only n^{k) (orig- 
inating from /5(fc) (dH)), but also n2{k) (related to f2{k) 
([251) 1 has to be included partially in the particle-state 
contribution n^z(N)^^^ momentum distribution. 

The latter leads to a successful explanation of the experi- 
mental data for the quasielastic scaling function f^^{'4>') 
(see Fig. 4 of Ref. 63]) extracted from inclusive electron 
scattering off nuclei. 

Finally, the normalized to unity proton (neutron) mo- 
mentum distribution [Eqs. ([T5)) - ([2^ ] has the form: 
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nz{N){k) = 



C{k)nziN){k) + Z{N)p[n2{k) + n^jk)] 
J dfc'{c(fc')^zw(fc') + Z{N)p[n2{k') + n^ik')] } 



(30) 



r 



(with normalization of nz(N)ik) presented by Eq. (US])). 

For the value of the parameter /3 we choose /3 ~ 5.0 
because of three reasons. Firstly, our calculations of 
nz(N)ik) for various nuclei using this value lead to re- 
sults very close to the empirical data for the nucleon mo- 
mentum distribution new (see Fig. 2 of Ref. [6l| and 
Fig. 3 of Ref. [63]) extracted in Refs. [H, [13 from the 
y-scaling analyses of inclusive electron scattering off nu- 
clei. Secondly, our estimations of the high-momentum 
components oi nz(N){k) [Eq. ((30|) ] showed that the value 
of (3 must be similar to the value of [Eq. ([^ J. It 
is shown (see Table I in [s^) that the value of es- 
timated from the height of the plateau exhibited by the 
ratio of the nucleon momentum distribution of a nucleus 
to the one of the deuteron at k > 2.0 fm~^ ranges from 
4.0 for to 4.4 for "OCa, 4.5 for ^^Fe, 4.8 for ^ospb, 
and 4.9 for nuclear matter {A=oo). Third, our results 
for the NMD (e.g. for ^^C, ^"^Wi and others) in LED ap- 
proach for large values of fc (fc > 2.0 fm^^) are similar to 
those obtained in the local density approximation using 
the Jastrow correlation method [80] , 



2. Theoretical approach based on the local density 
approximation 

It is well known that the inclusion of correlations in 
nuclear matter modifies the occupation probability pre- 
dicted by the local density Fermi gas model. It was shown 
in (soj that for a finite nucleus the separation between 
the mean-field contribution and correlation effects can be 
performed in an analogous way. According to Ref. [131 
one can introduce proton (neutron) momentum distribu- 
tion in a general way: 

nziN){k) - nfw (fc) + Snz^N){k), (31) 

where n^z(N)^^) mean- field contribution which can 

be momentum distribution corresponding to a Slater de- 
terminant generated by the single-particle wave functions 
in the momentum space or HE momentum distribution, 
while 5nz{N){k) embodies the corrections due to dynam- 
ical correlations not included in the MFA. If one applies 
the LDA to the second term of Eq. (|31|). the nucleon 
momentum distribution nz(N){k) can be written in the 
form: 

1 



rizw(fc)=nf(iv)(fc) 



47r3 



5v{kl^^\r),k)df, (32) 



where Si'{kp^'^\r), k) corresponds to the occupation 
probability that is entirely due to the effects of dynami- 
cal correlations induced by the NN interaction. The local 



Fermi momentum fc^*-^^ (r) is related to the proton (neu- 
tron) density through the relation 



fc^W(r) = [3n'pziN)ir)] 



1/3 



By definition of kp^^\r) one has / 5v{k^^'^\r), k)dk 
0. For convenience, a phenomenological procedure based 
on the results of the lowest order cluster (LOG) approx- 
imation developed in 'ss] has been followed to evaluate 
explicitly the correlated term. Choosing a correlation 
function of the form 

/(r) = 1 - e-^''^^ (34) 

the LOG gives for 5v{kP^\k) ^ 

5u{kl^''\T),k) - [Y{k,%)~kd,r]Q{kl^''\r)-k) 

+ 8{fcd„.r(fc,2)- [y(fc,4)]2}, (35) 

where 



(33) 



e + — e - 



with 



2k 



sgn(A;_) 



1 //.tNj3/2 



e 
/o 



dy (36) 



BVtt V2 



k± 



kp zL k X 



(37) 



(38) 



The quantity 



2[4^^f 
37r2 

1 / k 



[/(r)-l]V 



Z{N)' 



3V27r V 7 



(39) 



is the direct part of the Jastrow wound parameter. 

As in the approach based on the LED method dis- 
cussed before, in our work we use for n^^^'^(fc) the mo- 
mentum distributions obtained from DDHF calculations 



.DDHF 
''Z(N) 



(k) with normalization given by: 



n^l^^yik)dk^ZiN). 



(40) 
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For the densities pz{N){f) entering Eg. (p3)) we use the 
HF+BCS proton (neutron) densities [531 normahzed as: 

/ pf(X(r)dr = ^(iV). (41) 

For the correlation factor in Eq. p4|) we adopt the same 
value 7=1.1 fm^^ as in Ref. [80] which is taken from the 
microscopic nuclear matter calculations [s^ but seems 
to agree also with the data on momentum distributions 
n{k) in finite nuclei. 

III. RESULTS OF CALCULATIONS AND 
DISCUSSION 

We start by showing the results of our calculations 
for the monopole components of n(fc) and p(f) in the 
ioo,i2oa36gjj isotopes. As discussed in Ref. [H, [111 these 
are the only important components for the HF momen- 
tum distributions {n{k) ^ no(fc) = n{k) [Eq. p6|) ]V For 
easier reading, in the figures we will omit the subscript 
Z(A'^)-indices of n{k) for the proton and neutron momen- 
tum distributions replacing them by p{n), respectively. 
In Fig.[ljL are given the total momentum distributions for 
^"^Sn, ^^°Sn and ^^^Sn plotted in non-logarithmic scale. 
One can observe in the figure an appreciable difference 
between the curves in the low momentum region and, 
at the same time on this scale, no sensitivity of n{k) at 
fc > 1 fm^^ when increasing the number of neutrons. As 
can be expected, the momentum distributions calculated 
within the DDHF approach show a steep decrease up to 
fc < 2 fm-i that can be seen in all mean-field calculations. 
Here we emphasize that the effects of NN correlations on 
the momentum distributions can be seen at higher mo- 
menta (fc > 2 fm~^) when we will compare the DDHF 
results with those obtained in the approaches which take 
them into account. The results for n(fc) are related with 
the total density profiles of the selected three Sn isotopes 
shown in Fig. [T|d. For a comparison, in the same figure 
we present the total density p{r) of ^"°Sn obtained in 
Ref. 90] within the Hartree-Fock-Bogoliubov model us- 
ing the SLy4 parameterization of the Skyrme force. We 
note, following the discussion concerning Fig. 5 in [ssj . 
that the effect of adding more and more neutrons leads 
to an extension of the total density p{r) as one goes from 
^'^'^Sn to ^"^^Sn and to an emergence of a region at the 
nuclear surface quantified as a "neutron skin". The lat- 
ter is due to the larger spatial extension of the neutron 
density relatively to the proton one. The same consistent 
trends of the matter densities when increasing the mass 
number A were found in [o^] • One can see from Fig. [1)3 
that the densities of ^'"^Sn obtained in our previous work 
[53| and in Ref. [93] are very similar. Second, the same 
choice of two extreme neutron-deficient and neutron-rich 
isotopes and one stable isotope between them in the Sn 
chain 53| makes the difference in the results for n(fc) at 
low k more pronounced with respect to those for p{r) at 
large r. 



The effect of pairing correlations on the momentum 
distribution can be seen in Fig. [21 We restrict the discus- 
sion to the stable ^''Kr isotope which is weakly deformed. 
In this case the pairing correlations have been included 
by solving the BCS equations in the constant pairing gap 
mode. Fig.[2|shows that the effect of pairing correlations 
is small. At high fc this effect plays some role result- 
ing in more pronounced tails when BCS-correlations are 
included in the calculations. Nevertheless, as we will 
illustrate and discuss later, the effects of pairing cor- 
relations on the HF momentum distribution of nuclear 
matter and of finite nuclei are very different being much 
larger in the case of NM. Since there is no big differ- 
ence between the results for np(fc), n„(fc), and n(fc) with 
or without BCS-correlations included, afterwards we will 
use in our consideration only the results from the com- 
plete DDHF+BCS calculations. 

Our next step is to present and discuss the results for 
the NMD's of exotic nuclei obtained also in correlation 
methods. A comparison of these results for the neutron 
and proton momentum distributions of ^"^Ni, ^'^Kr, and 
^^°Sn nuclei is given in Fig. [3] together with the HF mo- 
mentum distributions. As can be seen, for all nuclei the 
inclusion of NN correlations strongly affects the high- 
momentum region of NMD. At fc > 1.5 fm^^ both LED 
and LDA momentum distributions start to deviate from 
the DDHF-I-BCS case. They behave rather similar in 
the interval 1.5 < fc < 3 fm~^. At fc > 3 fm~^ the 
LED method predicts systematically higher momentum 
components compared to LDA momentum distributions. 
This observation can be explained by the different ex- 
tent to which NN correlations are taken into account 
in both approaches. Our results for the NMD's in the 
LED method for large values of fc (fc > 2 fm^-'^) are 
similar to those obtained within the Jastrow correlation 
method and, thus, the high-momentum tails of n{k) are 
caused by the short-range NN correlations. The LDA ap- 
proach through the nuclear matter dynamic effects and 
using the local Fermi momentum k^^'^^ (r) calculated 
self-consistently by means of the HF density [Eq. (|33p ] 
produces less pronounced high-momentum tail, but still 
the results are very close to those obtained in the LED 
method. As was already shown, at fc > 1.5 fm~^ the 
DDHF-I-HF momentum distributions fall off rapidly by 
several orders of magnitude in contrast to the correlated 
NMD's. In addition, we observe that: i) the results 
shown above are similar for all nuclei in a given iso- 
topic chain and going from Ni to Sn isotopes, as well; 
ii) the behavior of n{k) is similar for protons and neu- 
trons; iii) at high fc the proton and neutron NMD's ob- 
tained within the LED method cannot be distinguished 
from each other because the high-momentum tails in this 
approach are determined by the high-momentum compo- 
nent of the nucleons in the deuteron [78| : iv) concerning 
the NMD's calculated in the LDA approach, some dif- 
ference between n(fc) for protons and neutrons can be 
observed due to Z( A^)-dependence of the local Fermi mo- 
mentum kp- 
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FIG. 1: (a) Total momentum distribution for ^""Sn, ^'^"Sn, and ^'"'Sn isotopes resulting from DDHF calculations. The nor- 
malization is: J n{k)dk = 1; (b) Total density distributions for the same Sn isotopes. The normalization is: / p{r)df = A. 
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In the next figures |4][7] we show the neutron n„(/c) 
[Figs, m and [B] and proton np(fc) [Figs. [5] and [7] mo- 
mentum distributions of some selected isotopes in the 
Ni and Sn chains, the same which have been consid- 
ered in Ref. [53| to calculate important nuclear properties 
in coordinate space. The results are presented in both 
logarithmic and linear scales in order to study the iso- 
topic sensitivity of these momentum distributions in the 
high-momentum region and in the region of small mo- 
menta, respectively. In addition, in each of the figures 
the results for neutron and proton momentum distribu- 
tions in the DDHF-I-BCS method and in the correlation 
LFD and LDA approaches at fc < 2 fm~^ are given sep- 
arately in panels (b) and (c). First, it is seen from the 
figures that the general trend in the behavior of NMD's 
obtained within the methods used in the calculations and 
already shown in Fig. [3] is preserved. Second, the evolu- 
tion of the NMD's as we increase the number of neutrons 
consists of an increase of the high-momentum tails (for 
k > 1.5 fm~^) of nn{k), while the effect on np{k) is op- 
posite. However, the spreading of the tails corresponding 
to np{k) of the considered isotopes is of the same order 
although the number of protons remains the same. In 
this respect, the results presented in Figs. O and [7] are 
challenging because they show how proton momentum 
distributions "feel" the different number of neutrons in 
exotic nuclei. We would like also to emphasize that the 
LFD method does not show this isotopic sensitivity, in 
contrast to the HF and LDA methods which still demon- 
strate this trend. Concerning the low-momentum region 
it can be seen from Figs. |4][7] that NMD's are very sensi- 
tive to the details of the calculations. In this region n„(fc) 
decreases while, on the contrary, np{k) increases with the 
increase of the number of neutrons N. This is a common 
feature of the calculated results obtained in all methods. 
Nevertheless, in this region the spreading is considerably 
reduced. From the comparison of the proton and neu- 



tron momentum distributions at low momenta it comes 
out that the protons exhibit the same trend [Figs. [5] and 
[71[(b) and (c))], while the neutrons are more sensitive to 
nuclear shell effects (see panels (b) and (c) of Figs.|4]and 
HI). 

The isotonic sensitivity of the neutron, proton and to- 
tal momentum distributions of "^^Ni, *^Kr, and ^'^°Sn nu- 
clei (A^=50) is shown in Fig. [51 It can be seen from 
the figure a small difference between the different curves 
when using a given approach (DDHF-fBCS or LDA). 
The neutron momentum distributions n„(fc) do not differ 
so much in comparison to the proton momentum distri- 
butions np{k). The relative contributions of n„(fc) and 
np{k) to the total momentum distribution n{k) lead to 
almost equal high-momentum tails of n(fc) for these iso- 
tones. This is in accordance with the well-known results 
showing that the tails of n{k)/A are almost equal for all 
nuclei (for a review, see e.g. [64] ) . 

The role of the deformation of the mean field on the 
NMD is studied on the example of ^^Kr isotope. For 
this purpose, we show in Fig. [9] the intrinsic momentum 
distribution for neutrons and protons and for oblate and 
prolate shapes of ^^Kr. The results of the three methods 
considered in our work are given and compared together 
with the result for the spherical case obtained within the 
HF method. The differences between the momentum dis- 
tributions calculated within a given theoretical method 
are negligible (especially at high k) and practically can 
not be distinguished. Thus, we find almost no depen- 
dence of the ri„(fc) and np{k) on the character of defor- 
mation. 

The effects of the correlations included in the LFD and 
LDA methods on the neutron momentum distribution of 
^^°Sn are presented in Fig. [TOb . The parameters that 
govern the correlations in the two approaches are (3 and 
7, respectively. The results of the calculations for three 
values of each of them are shown by thick lines in the case 
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FIG. 4: (a) Neutron momentum distributions obtained within the DDHF+BCS, LED, and LDA methods for ^"Ni (solid line), 
^*Ni (dashed line), and ^*Ni (dotted line) isotopes. The normalization is: J nn{k)dk = 1. The DDHF+BCS results, as well as 
the LFD and LDA results are separately shown in a linear scale in (b) and (c), respectively. 
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of LDA and by thin lines in the case of LFD method. The 
value (3=5 is the upper limit of this parameter [s^ and, 
therefore, we give in Fig. [TUk also the results for n„(fc) for 
two smaller values /3=4.4 and /3=4.7. It turns out that 
ri„(fc) does not change significantly for different values 
of /3, thus showing the strong presence of correlations at 



short distances within the LFD method not only for the 
stable, but also for the exotic nuclei. However, a larger 
sensitivity of ri„(fc) on the parameter 7 in the LDA ap- 
proach appears , particularly in the interval 1.5 < A: < 3 
fm~^. The value 7=1.1 fm^^ (from Rcf. [80]) and two 
more values 7=0.9 fm^^ and 7=1.3 fm^^ are used in the 
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calculations. As expected, the obtained neutron momen- 
tum distributions contain smaller correlation effects at 
larger values of 7. This is in accordance with the behav- 
ior of the correlation function /(r) (see Eq. ([M]) ). At 
fc > 3 fm^^ the LDA results start to deviate in an op- 
posite way, but in this very high-momentum region no 
definite conclusion about the role of correlations can be 
drawn. 

In our work we study also the sensitivity of our results 
for NMD to different effective NN forces. In Fig. [TUb we 
show the neutron momentum distributions of -'^^"Sn ob- 
tained by using SLy4 Skyrme force together with the re- 
sults obtained from other parametrizations, namely, Sk3 
[qH and SG2 [o^] . For stable spherical nuclei it is known 
that all the Skyrme-type effective interactions give simi- 
lar results for the total momentum distribution ,93J. We 
would like to resolve possible ambiguities concerning ex- 
otic nuclei by comparing the momentum distributions for 
neutrons on the example of the stable ^^"Sn isotope. In 
fact, some sensitivity to different Skyrme forces can be 
observed in the inset of Fig. [TUb for momenta k < 0.5 
fm~^. This is due to shell effects that are found also 
in the density profiles of the nucleus. Generally, as it 
has been pointed out in Ref. [s^, the predictions for the 
charge and matter rms radii do not differ too much when 
different Skyrme forces are used. In analogy, different 
Skyrme interactions also produce similar results for the 
NMD. 

Finally, in Fig. [TT] a comparison of the results for the 
total momentum distribution n{k) of nuclear matter cal- 
culated within the MFA and correlation methods used 
in our work is shown. The HF momentum distribution 
is strongly affected by pairing correlations which build 
up a long tail at high momentum (fc > kp). Comparing 
this result with the result illustrated in Fig. [2] for ^^Kr, 
the different role played by pairing correlations on the 
DDHF momentum distribution of NM and of finite nu- 
clei becomes clear. Moreover, it is interesting to explore 
the case when one includes other type of correlations in 
NM. Stringari et al. [s^l have already shown in their 
model based on the LDA the prediction for n{k) in the 
case of nuclear matter. An enhancement of the high- 
momentum components of n(k) can be seen from Fig. 1111 
when both LDA and LFD methods are used. Hence, in 
nuclear matter the effects of short-range and tensor cor- 
relations are much stronger than the BCS-correlations 
taken into account in the DDHF-I-BCS calculations. 



IV. CONCLUSIONS 

In this work we investigated the properties of the mo- 
mentum distributions of medium and heavy exotic nu- 
clei, especially of Ni, Kr, and Sn even-even isotopes. The 
theoretical study was performed on the base of the mean- 
field method, as well as of two correlation methods tak- 
ing into account the NN correlations at short distances. 
The neutron, proton, and total momentum distributions 



of these nuclei were calculated within: i) a deformed 
DDHF-I-BCS method with Skyrme-type effective inter- 
actions [13, nil; ii) a theoretical approac h |6ll . [63l pTSj 
based on the light-front dynamics method [79] ; iii) a the- 
oretical model based on the local density approximation 
[sot . I'^ the DDHF+BCS calculations we consider the 
monopole component of n{k) {n{k) = no{k)) since this 
is the only important component in the expansion of the 
HF ground-state momentum distribution [Eq. p6p ]. On 
the other hand, the two correlation approaches allow one 
to include both the mean-field and short-range effects for 
the description of the nucleon momentum distribution. 

The study of the isotopic sensitivity of various kinds 
of momentum distributions shows different trends. For a 
given isotopic chain, we find that in the high-momentum 
region (k > 1.5 fm~^) the high- momentum tails of the 
neutron momentum distributions nn{k) increase with the 
increase of the number of neutrons N, while the proton 
momentum distributions np{k) exhibit an opposite ef- 
fect. In the same region the LFD method does not show 
this isotopic sensitivity, in contrast to the DDHF-I-BCS 
and LDA methods. At low momenta n„(fc) decreases 
while, on the contrary, np{k) increases with the increase 
of N . Additionally to the isotopic sensitivity we studied 
how the momentum distributions of some isotones are 
modified keeping the neutron number constant. We find 
that the total momentum distributions of ^^Ni, ^^Kr, and 
^°°Sn nuclei (iV=50) reveal the same high- momentum 
tails in all methods used. 

The role of the deformation on the momentum distri- 
butions is discussed in the present work on the example 
of ^^Kr isotope. As it has been found in Refs. [4l,[93|, the 
isotropy of the total momentum distribution is a prop- 
erty of the nucleus at equilibrium. Our results for the 
neutron and proton momentum distributions of this nu- 
cleus show small changes in the overall behavior for the 
oblate and prolate shapes. Although the neutron and 
proton densities change with deformation [ssj , the mo- 
mentum distributions demonstrate a very weak depen- 
dence on the character of deformation. This is valid for 
all three theoretical approaches explored in our work. 

The pairing correlations are shown to influence the 
high- momentum behavior of the neutron n„(fc), proton 
rip(fc), and total n(k) momentum distributions in the case 
of *'*Kr, but the differences between the results with or 
without BCS-correlations included in the calculations are 
very small. The latter are even smaller considering nuclei 
from the Ni and Sn isotopic chains that have a spheri- 
cal shape. The effect of pairing correlations on the HF 
momentum distribution is much stronger in the case of 
nuclear matter producing a tail for momenta k > kp. 

The effects of the dynamical correlations on the mo- 
mentum distributions have been calculated using ap- 
proaches which account for correlations at short dis- 
tances. In the correlation approach based on LFD, the 
incorporation of the LFD result for the nucleon momen- 
tum distribution in the deuteron makes it possible to 
study NMD's also for exotic nuclei, especially at high 
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FIG. 10: (a) Neutron momentum distributions obtained within the LFD (for different values of the parameter /3) and LDA (for 
different values of the parameter 7) methods for ^^'^Sn; (b) Neutron momentum distributions obtained within the DDHF+BCS 
method for ^^"Sn and for different Skyrme forces. The normalization is: / n„(fc)dfc = 1. 
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FIG. 11: Comparison of DDHF results for the momentum distributions of nuclear matter with (dashed line) and without (solid 
line) pairing correlations with the results from LFD (dotted line) and LDA (dash-dotted line) methods. 
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momenta (A; > 2 fm~^). It is known [H, that the 
latter are rescaled versions of n(fc) in the deuteron. This 
is the main reason for the fact that the LFD calculations 
do not show isotopic sensitivity on the obtained NMD's. 
Concerning the NMD's calculated in the LDA approach 
in which a Jastrow-type correlation function is adopted, 
some difference between n{k) for protons and neutrons 
is observed due to 2'(A^)-dependence of the local Fermi 
momentum kp- In general, a strong enhancement of the 
momentum distributions at large k is observed in com- 
parison to the MFA result and, at the same time, both 
LFD and LDA results are similar in the high-momentum 
region. In our work we pay a particular attention to the 
dependence of the results for the momentum distribu- 
tions on the correlation parameters (3 and 7. In the LFD 
calculations we used the value of the parameter /3 = 5.0. 
It turned out that this value leads to a good agreement 
with the empirical data for n{k) extracted from the y- 
scaling analyses of inclusive electron scattering off nuclei 
(55I \8T\ . The LDA calculations have been carried out 
for 7=1.1 fm~^, the same value being adopted in nuclear 
matter and also providing a good choice for finite nuclei. 
In addition, we note that our prediction for n(k) in NM 
obtained within the LFD method is in agreement with 
the LDA result [13] • In our opinion, however, the ques- 
tion for the specific values of the parameters (3 and 7 
which determine the strength of the correlations is still 



open. 

We would like to emphasize that, in our work, a possi- 
ble practical way to make predictions for the momentum 
distributions of exotic nuclei far from the stability line is 
proposed that provides a systematic description of n{k) 
in medium- weight and heavy nuclei. The comparison of 
the predicted nucleon momentum distributions with the 
results of possible experiments using a colliding electron- 
exotic nuclei storage rings would show the effect of the 
neutron excess in these nuclei and will be also a test of 
various theoretical models of the structure of exotic nu- 
clei. 
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